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Mechanical oscillations of synchronous machines following the 
application of sudden shaft loads are important to operating engi- 
neers. This paper presents an analytical technique to determine the 
stability of a synchronous machine. We first decompose the nonlinear 
system representing the machine into an infinite number of subsys- 
tems connected in parallel. We then determine a bound on the input 
to the machine for which the number of subsystems will effectively be 
finite in number, then determine stability of the entire system by 
determining the stability of each subsystem. 

I. INTRODUCTION 

The problems associated with the maintenance of stability of syn- 
chronous machines and the factors affecting stability during transient 
disturbances have received considerable attention in recent years. If 
the sudden load application is such that synchronizing torque is less 
than the load torque, the machine will be thrown out of synchronism 
and instability will occur. A synchronous machine has steady-state 
stability if, after a small slow disturbance, it can regain and maintain 
synchronous speed. On the other hand, if the machine can regain and 
maintain synchronous speed after large sudden disturbances, it has 
transient stability. The maximum steady-state operating limit will 
depend upon the magnitude of instantaneous load changes and its 
ability to follow quickly any load changes. It is the purpose of this 
paper to determine an upper bound on the input load so that synchron- 
ism will never be lost. 

II. MATHEMATICAL MODEL 

The equation of motion of the synchronous machine is well known 
and has been extensively discussed in the literature. 1-6 However, the 
derived equation based on the electromechanical properties of the 
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system is highly nonlinear, and the determination of the time-domain 
response of the system seems to be prohibitively difficult. Thus, to 
represent the system by some suitable mathematical model, the follow- 
ing simplifying assumptions are made. 

(i) The twist in the shaft is negligibly small since the length of the 
shaft is small and the rotational stiffness is large. 

(ii) All inductances are independent of current since nonsaturated 
conditions are considered. 

(Hi) The armature flux wave is sinusoidally distributed. 

(iv) The currents in both the armature and field circuits are consid- 
ered to be independent of a sudden change in operating conditions. 
Results obtained with an induction motor having a rotor winding 
excited with direct current appear to justify this assumption for those 
cases in which the time constants of the stator and rotor windings are 
relatively small or the moment of inertia is large enough to assume 
small acceleration. 

(v) All resistances are negligible compared to inductances. 

(iv) The direct-axis synchronous reactance is equal to the quadra- 
ture-axis synchronous reactance. This assumption is quite valid for a 
nonsalient pole machine but for a salient pole machine this assumption 
is only approximate. 

In order to write the system equations, the following symbols are 
introduced. 

J = polar moment of inertia of the rotor. 

w = instantaneous speed. 
<o s = synchronous speed. 

8 = angle between the rotor pole position on load 

and its positions on no load (rotor angle). 
Ti = externally applied torque. 
T e = torque of electrical origin. 
D = angular viscous friction coefficient. 
Pi = electrical power input. 
Pd = system damping factor. 
P m = maximum static synchronous power. 
P = output power. 
M = inertia constant. 

Based upon the simplifying assumptions, the equation of motion of 
the synchronous machine can now be derived. However, the detailed 
derivation will not be given here, since it can be found elsewhere. 1-7 
Only an outline of the derivation which is necessary for the subsequent 
discussion is presented below. 

8 = (a - us)t, (1) 
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— = 03 - 03 s , \£) 

at 
d 2 8 do3 , Q . 

~d? = Tf (3) 

The equation of motion is given in Ref. 7 as 

Ti = J^ + Dto + Te. (4) 

ar 

By substituting (2) and (3) in (4) and by converting the torque equation 
(4) into power equation 

d 2 8 _ (d8 



or 



ft— /y+n^ s +*l + ''. 



where P d (dS/dt) = Du(d8/dt) + Do303 s , Pd is the system damping factor 
and can only be determined empirically. P a has been evaluated 6 as P m 
sin 8. Equation (5) can now be put into dimensionless form by intro- 
ducing the following change in variables: 

Pi 

P = TT 



2£ = 



J{MP m ) 
By using the above transformation, (5) becomes 



<?8 0> d8 
or 



- ri + 2£ -j- + sin 8 = p 



8 + 2£8 + 8 + N(8) = p, (6) 

where N(8) is a nonlinear function of 8 and is given by 

The behavior of the synchronous machine can now be determined 
from (6). The damping coefficient £ cannot be evaluated using exact 
methods. Its value is usually approximated empirically from the knowl- 
edge of the system behavior. A value of < £ < 0.4 seems to be a 
reasonable choice. 
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The input-output relationship of the synchronous machine can now 
be viewed in the form of a simplified block diagram as shown in Fig. 
1. The input load ratio p(r) is the input to the system, and its response 
is the rotor angle 5(t). Now by functional expansion, 8 8(t) can be 
expressed in terms of p(r) with all the limits of integrations ranging 
from to oo 

5(t) = I hAtJpir - h) dti + I h 2 (t h t 2 )p(T - ti)p(T - t 2 ) dti dt 2 
+ 1(1 h 3 (t u t 2 , t 3 )p(r - h)p(r - t 2 )p(T - t 3 ) dti dt 2 dt 3 + • • • . (7) 



Let 

5„(t) =1 • • • I hn(t h t 2 , • • • , t n )p(r - U) 



• • • p(r - t n ) dti • • • dt n . (8) 



Thus 



5(t) = I S„(t). 

n=l 



Hence the nonlinear system in Fig. 1 can be broken down into an 
infinite number of subsystems connected in parallel as shown in Fig. 
2. The response of the first-order or linear subsystem is 5i(t), and its 
impulse response is hi(ti). The response of the nth order subystem is 
5„(t), and its impulse response is h n (t\, t 2 , • • • , t n ). 

Thus the system shown in Fig. 1 will be stable if 

(i) Equation (7) is convergent, which also implies | 5, | > | 5,-+i | . 

(ii) The linear subsystem is stable, which also implies s 2 + 2£s + 1 
= does not have any root with positive real part, s = d/dt. 

Condition (ii) can easily be checked from the system equation. To 
prove the condition of convergence, we proceed as follows. 

Let h(t) be the impulse response of the first order or the linear 
subsystem, i.e., 

h(t) ~ 



and 



sr + 2 £s + 1 



H= \h(t)\dt. 



P(T] 



SYNCHRONOUS 
MACHINE 



8(T) 



Fig. 1— Input-output relationship of a synchronous machine. 
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<y, in 



P(T) 




SlT) 



Let 



Fig. 2 — Decomposition of a synchronous machine into infinite subsystems. 



max |p(r)| < P. 



With the above assumptions, by the technique of series reversion 
method, 9 Barrett has shown that the functional expansion expressed 
in (7) will be convergent if the input has an upper bound given by 

1 



I I IT 

|p(t) I < P = — =— cosh" 



— sinn cosh 



H H H 

The integral H = J | h(t) \ dt can be shown 12 to have the value 



H = coth 



2^/T^f 



By knowing H, P can be calculated. The values of £, H, and P are 
shown in Table I. 

It can also be shown 12 that, for the system described by (6), /12U1, t 2 ), 
h.\(t\, • • • , U) are zeros. Hence 

5(t) = 5,(t) + Ur) + &(t), ••-. 

III. DETERMINATION OF THE STEADY-STATE RESPONSE 

A sudden change in the load on the synchronous machine is equiv- 
alent to a step change in the prime mover torque applied to the shaft. 
Hence, the system stability can thus be investigated on the basis of a 
step input. It is also clear from Fig. 2 that the steady-state response of 
the entire system is the summation of the steady-state responses of 
each of the subsystems. Let 

A„ = lim 8„(t) 
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Table I — Input as a function of damping 
coefficient 



Damping 




Upper Bound 


Coefficient 


Jo \h(t)\dt 


of Input 


£ 


H 


P 


0.01 


63.664 


0.002 


0.02 


31.835 


0.005 


0.03 


21.227 


0.010 


0.04 


15.924 


0.015 


0.05 


12.743 


0.021 


0.06 


10.623 


0.027 


0.07 


9.109 


0.034 


0.08 


7.974 


0.042 


0.09 


7.092 


0.050 


0.10 


6.387 


0.058 


0.11 


5.810 


0.067 


0.12 


5.330 


0.076 


0.13 


4.924 


0.085 


0.14 


4.576 


0.095 


0.15 


4.275 


0.105 


0.16 


4.012 


0.116 


0.17 


3.780 


0.127 


0.18 


3.574 


0.138 


0.19 


3.390 


0.149 


0.20 


3.225 


0.160 


0.25 


2.599 


0.221 


0.30 


2.186 


0.285 


0.35 


1.895 


0.353 


0.40 


1.680 


0.421 


0.45 


1.517 


0.490 


0.50 


1.390 


0.557 



and the input to the system be a step function p(r) = pi. It has already 
been determined from Table I that, for stability the values p, will be in 
the range p, < P < 0.421 for < £ < 0.4. The values of Ai, A 3 , • • • can 
now be determined by the application of the final- value theorem. Ai is 
given as 

spi 



Ai = lim 



o s{s* + 2£s + 1) 



= Pi 



To determine An, we proceed as follows: First we take a three- 
dimensional transform 10 of the third subsystem corresponding to the 
output 5 3 (t). Next, we apply the technique of association of variables 11 
and make a transformation from three dimensions to a single dimen- 
sion. Now we are ready to apply the final- value theorem to find A3 . 

ha(ti, £2, t 3 ) <-> H 3 (s\, S2, s 3 ). 

Let A stand for associate variables 81, S2, s 3 — * s. Then 

H 3 (si, s 2 ,s 3 )p'f 



S 3 (s) = A 



S1S2S3 



A3 = lim sA 



H 3 (si, s 2 , s 3 )pf 



S1S2S3 
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This limit has been evaluated 12 to be 

b 

Higher limits A 5 , A7, etc., will further be smaller, and the evaluation of 
these terms is not necessary since their contributions will be insignifi- 
cant. Thus 

A = lim 5(t) = Ai + A3 

is a very good approximation, and 

max A = pi ( 1 + y J < v/2, 

since < p, < 0.421. Physically, it means that if the magnitude of the 
abrupt load change is such that < Pi/P m < 0.421, the system will 
never lose synchronism and will remain stable. Thus, for p(r) < 0.421 
the system can be described by the linear transfer function within a 
reasonable degree of accuracy. 

With the above assumption, the output is given by Ref. 13: 



8(r)= Pi 



1 - e~* T sin(Vl - £ 2 t + cos -1 © 



S( t) can b e expressed in terms of t only by using the relationship t = 
ry/M/K. 

Now by separating p into two parts, the initial load ratio p 00 and the 
abrupt load ratio pl , the maximum value of pl can be determined 
before the system goes out of synchronism. 
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